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Abstract 

We present predictions on the total cross sections and on the ratio of the real part to the 
imaginary part of the elastic amplitude (p parameter) for present and future pp andpp colliders, and 
on total cross sections for jp — > hadrons at cosmic-ray energies and for 77 — > hadrons up to y/s = 1 
TeV. These predictions are based on an extensive study of possible analytic parametrisations 
invoking the biggest hadronic dataset available at t = 0. The uncertainties on total cross sections, 
including the systematic errors due to contradictory data points from FNAL, can reach 1.9% 
at RHIC, 3.1% at the Tevatron, and 4.8% at the LHC, whereas those on the p parameter are 
respectively 5.4%, 5.2%, and 5.4%. 

PACS numbers: 13.85.-t, 11.55.-m, 12.40.Nn, 13.60.Hb 
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In recent works || , we have performed an exhaustive study of the analytic parametri- 
sations of soft data at t — 0. For this purpose, we gathered the largest available set of data 
at t — 0, which includes all measured total cross sections and ratios of the real part to the 
imaginary part of the elastic amplitude (p parameter) for the scattering of pp, pp, ^p, K ± p, 
and total cross sections for jp, 77 and £~p || |]. 

Several experiments are under way [f|], or being planned, to measure the hadronic am- 
plitudes at t — 0. Some authors || § also presented what they feel are reference values for 
the total 7p and 77 — > hadrons cross sections. Thus it is timely and appropriate to present 
independently our predictions for the forward observables at RHIC, the Tevatron-run II and 
the LHC as well as for 7p total cross section at cosmic-ray energies and for 77 total cross 
sections up to 1 TeV. 

We can summarize the general form of the parametrisations by quoting the form of total 
cross sections, from which the p parameter is obtained via analyticity. The ingredients are 
the contribution Y ab of the highest meson trajectories (p, u, a and /) and the rising C = +1 
term H ab from the pomeron contribution to the total cross section, which can be written for 
the scattering of a on b: 

o? ot = {Y ab + H ab )/s (1) 

The first term is parametrised via Regge theory, and we allow the lower trajectories to be 
partially degenerate, i.e. our experience shows that it is enough to introduce one intercept 
for the C — +1 trajectories, and another one for the C = — 1 0. A further lifting of the 
degeneracy is certainly possible, but does not seem to modify significantly the results ||. 
Hence we use 

yab = yab ^ J ^ ± Y f ( S / (2) 

with Si = 1 GeV 2 . The contribution of these trajectories is represented by RR in the model 
abbreviations. 

As for the part rising with energy, we consider here two main options: it can rise as a 
logs, or as a log 2 s, with in each case the possibility to add a constant term. We shall not 
consider the simple-pole parametrisation , not only because it is disfavored by our ranking 
procedure (see below), but also because we want to make predictions at very high energies, 
where unitarisation must set in (JB]j . In the following, we shall only refer explicitly to our 
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preferred parametrisation of H ab , which we note as PL2: 

H ab = s{B ab In 2 {s/ s ) + P ab ) (3) 

where so is a universal scale parameter (to be determined by the fits) identical for all colli- 
sions. 

We have considered several possible constraints on the parameters of Eqs. (@-|3]): degen- 
eracy of the reggeon trajectories (a+ = a_); universality of rising terms (B ab independent 



of the hadrons) |10|, |TT|; factorization for the residues in the case of the 77 and 7p cross 
sections (i? 77 = SH 1P = 5 2 H pp ); quark counting rules |L2| (predicting the cross section 
from pp, Kp and 7rp); and finally the Johnson- Treiman-Freund [DJ relation for the cross 
section differences. 

Out of the 256 possible variants, we showed that 24 met our criteria for applicability (an 
overall x 2 /dof < 1.0 and a non-negative pomeron contribution at all energies) if one fitted 
only to a tot for y/s > 10 GeV, and 5 did for y/s > 4 GeV (see Table XI from []]]), whereas 20 
(resp. 4) variants obeyed this criterion when a fit to both <7 to t an d p was performed, for y/s > 
10 (resp. 5) GeV (see Table XIV from JE). We shall neither give here the list of models, nor 
spell out ranking criteria based on new indicators that quantify certain qualities of the fits, 
but simply mention that the triple-pole parametrisation RRP n /L2 u [|l(], |Tl| was determined 
to be the highest-ranking model leading to the most satisfactory description of the data (see 
similar conclusions in Jl4|). This parameterization has a universal (u) B\og 2 (s/so) term, a 
non-factorizing (nf) constant term and non-degenerate lower trajectories. 

We start by giving the predictions of this model, adjusted for (y/s > 5 GeV), with 
updated data points from ZEUS |L5] . These predictions include statistical errors calculated 
from the full error matrix Ey. We define 

with Q = o to t or p and Xi the parameters of the model. These errors are shown in Figs. 1 
and 2 by a filled band, and in Tables fl], |T|, and [IV. 



In these figures and tables, we also give our estimate of the systematic uncertainty coming 
from the discrepancy between different FNAL measurements of a tot- we fit RRP n jL2 u either 
to the high data (CDF) or to the low ones (E710/E811), and get two error bands. The 
distances from the central value of the combined fit to the upper (resp. lower) border of 
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FIG. 1: Predictions for total cross sections. The black error band shows the statistical errors to 
the best fit, the closest curves near it give the sum of statistical and systematic errors to the best 
fit due to the ambiguity in Tevatron data , and the highest and lowest curves show the total errors 
bands from all models considered in this letter (note that the upper curve showing the systematic 
error is indistinguishable from the highest curve in this case). 
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FIG. 2: Predictions for the p parameter. The curves and band are as in Fig. 1. 

these bands give us the positive (resp. negative) systematic errors. We estimate the total 
errors as the sum of the systematic and of the statistical uncertainties . 

One can see that the total errors on total cross sections are of the order of 1.9% at RHIC, 
of the order of 3.1% at the Tevatron and as large as 4.8% at the LHC and dominated by the 
systematic errors. The errors on the p parameter are much larger, reaching 5.4% at RHIC, 
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5.2% at the Tevatron and 5.4% at the LHC This is due to the fact that experimental errors 
are bigger, hence less constraining, but this also stems from the incompatibility of some 
low-energy determinations of p 0. This means that the systematic error is always bigger 
than the statistical one. 

Concerning the contradictory data, we are forced to use them in our fits until the discrep- 
ancy is resolved by further experiments. In the case of the Tevatron data, one can see that 
the discrepancy results in a big shift (of more than la) in the central value of the coefficient 
B of the log 2 s term, which controls the asymptotic behavior, and hence that asymptotic 
predictions are appreciably weakened by the present situation. The opportunities to measure 
a tot and p will be scarce in the future, hence any new measurement at RHIC, the Tevatron 
run II and the LHC should not be missed. Unfortunately, the recent publication of E-811 
T6fl does not clear the problem as their value for p is fully compatible with our preferred 



model, whereas their number for a tot (which is highly correlated with p) has hardly changed. 

It is interesting to note that the choice of one FNAL result or the others leads to a 
variation of the overall fit quality, as shown in Table | (last two columns) Q: the variant 
with CDF data has slightly better overall x 2 /dof and better x 2 jnop distribution over sub- 
samples. We can consider this as an indication that the global picture emerging from fits to 
all data on forward observables supports the CDF data and disfavors the E710/E811 data 
at y/s = 1.8 TeV (see also an analogous conclusion based on the other arguments in JT7[). 

Finally, we also present in Figs. 1 and 2 our estimate of the region where new physics 
would be discovered. For each of the 20 parametrisations which satisfy our criteria for 
applicability [Q for y/s > 10 GeV, and which obey the Froissart-Martin bound fl8], ||, we 
construct errors bands according to Eq. (Q). This gives us 20 lcr-error bands. Their union 
represents the "allowed region" where analytic models built according to (|l|) can reproduce 
the data. A measurement outside of this region would imply that new physics ingredients 
are needed. 

To conclude, we believe that we have given here the best possible estimates for present 
and future pp and pp facilities. Although one might be tempted to use only data in an 
energy range close to the one measured, one must realize that analytic parametrisations 
are constrained both by lower-energy data, and by their asymptotic regime. Because the 
pomeron mixes (physically and numerically) with the / trajectory, fits to all data help to 
disentangle the two contributions. 
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TABLE I: Summary of the quality of the fits at different stages of the Review of Particle Physics 
(RPP) database (DB): DB02 - The 2002 RPP DB; DB02Z - the 2002 RPP DB with new ZEUS 
data, DB02Z-CDF - with the CDF point removed; DB02Z-E710/E811 with E710/E811 points 
removed. The first line gives the overall x 2 /dof for the global fits, the other lines give the x 2 /nop 
for data sub-samples, the last line gives in each case the parameter controlling the asymptotic form 
of cross sections. 





DB02 


DB02Z 


DB02Z 


DB02Z 


Sample 






-CDF 


-E710/E811 


total 


0.968 


0.966 


0.964 


0.951 


total cross sections 


pp 


1.15 


1.15 


1.12 


1.05 


pp 


0.84 


0.84 


0.84 


0.84 


1T~p 


0.96 


0.96 


0.96 


0.96 


7T + p 


0.71 


0.71 


0.71 


0.71 


K-p 


0.62 


0.62 


0.62 


0.61 


K+p 


0.71 


0.71 


0.71 


0.71 


T,~p 


0.38 


0.38 


0.38 


0.38 


IP 


0.61 


0.58 


0.58 


0.58 


77 


0.65 


0.64 


0.64 


0.63 


elastic forward Re/Im 


pp 


0.52 


0.52 


0.52 


0.53 


pp 


1.83 


1.83 


1.83 


1.80 


1T~p 


1.10 


1.10 


1.09 


1.14 


17+P 


1.50 


1.50 


1.52 


1.46 


K-p 


1.00 


0.99 


1.01 


0.96 


K+p 


1.07 


1.07 


1.10 


0.98 



values of the parameter B 



0.307(10) 0.307(10) 0.301(10) 0.327(10) 
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A sharpening of our error bars would enable one to decide if the unitarisation plays an 
essential role and what form it takes. This in turn can have an impact on the determination 
of the survival of probability gaps in hard scattering, and on the usefulness of pomeron 
exchange clS db detection tool. 

Any significant deviation from the predictions based on model RRP n /L2 M will lead to a re- 
evaluation of the hierarchy of models and presumably change the preferred parametrisation 
to another one. A deviation from the "allowed region" would be an indication that strong 
interactions demand a generalization of the analytic models discussed so far, e.g. by adding 
odderon terms, or new pomeron terms, as suggested by QCD. 
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